Divergence of the Griineisen Ratio at Quantum Critical Points in 

Heavy Fermion Metals 

R. KiichlerW, N. Oeschler^, P. Ge genwart^ 1 -*, T. Cichorek^ 1 ), K. Neumaier^ 2 \ 
0. Tegus (3) , C. Geibel (1) , J.A. Mydosh«'( 4 \ F. SteglichW, L. Zhu (5) , and Q. Si (5) 
^Max-Planck Institute for Chemical Physics of Solids, D-01187 Dresden, Germany 
( 2 > Walther Meissner Institute, D-857^8 Garching, Germany 
® Van der Waals-Zeeman Laboratory, 
University of Amsterdam, The Netherlands 
W Kamerlingh Onnes Laboratory, Leiden University, The Netherlands 
^Department of Physics and Astronomy, 
Rice University, Houston, TX 77005-1892 

Abstract 

We present low-temperature volume thermal expansion, (3, and specific heat, C, measurements 
on high-quality single crystals of CeNi2Ge2 and YbRli2(Sio.95Geo.o5)2 which are located very near 
to quantum critical points. For both systems, (3 shows a more singular temperature dependence 
than C, and thus the Griineisen ratio r oc (3/C diverges as T — » 0. For CeNi2Ge2, our results 
are in accordance with the spin-density wave (SDW) scenario for three-dimensional critical spin- 
fluctuations. By contrast, the observed singularity in YbRli2(Sio.95Geo.o5)2 cannot be explained by 
the itinerant SDW theory but is qualitatively consistent with a locally quantum critical picture. 

PACS numbers: 71.10.HF,71.27.+a 
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Quantum critical points (QCPs) are of extensive current interest to the physics of corre- 
lated electrons, as proximity to a QCP provides a route towards non-Fermi liquid behavior. 
While a broad range of correlated electron materials are being studied in this context, heavy 
fermions have been playing an especially important role: a growing list of heavy fermion 
(HF) metals explicitly displays magnetic QCPs Systematic experiments in 

these systems promise to shed considerable light on the general physics of quantum critical 
metals. Indeed, recent experiments have shown that, at least in some of the HF metals, 
the traditional theory of metallic magnetic quantum phase transition fails. This traditional 
picture jf| describes a T = SDW transition and, relatedly, a mean-field type of quantum 
critical behavior. More recently, it has been shown that a destruction of Kondo resonances 
can lead to a breakdown of the SDW picture [?], Q ; what emerges instead are new classes of 
QCPs that are locally critical []]]. 

Given these experimental and theoretical developments, it seems timely to address the 
conditions under which these different types of QCPs arise. For this purpose, it would be im- 
portant to carry out comparative studies of different heavy fermion materials. This paper re- 
ports one such study. We have chosen the HF systems CeNi 2 Ge2 and YbRh 2 (Sio.95Ge .o5)2 
1^,0], both of which crystallize in the tetragonal ThCr 2 Si2 structure. Both are ideally suited 
to study antiferromagnetic (AF) QCPs since they are located very near to the magnetic 
instability, and since the effect of disorder is minimized in these high quality single crystals 
with low residual resistivities. We have focused on measurements of the thermal expansion, 
P, and Griineisen ratio, r oc P/C, where C denotes the specific heat, since recent theoretical 
work has shown that i) T is divergent as T goes to zero at any QCP and ii) the associated 
critical exponent can be used to differentiate between different types of QCP. 

Presently measurements of the thermal expansion and Griineisen ratio for systems located 
directly at the QCP are lacking. Only for Ce 1 _ a ,La a ;Ru2Si2, which orders antiferromagneti- 
cally for x > x c with x c = 0.075, P has been measured for concentrations x = and x = 0.05 
at temperatures above 0.4 K. A very large T was obtained which, however, was found to 
saturate at low temperatures [llj]. In other solids too, all previous measurements reported 
in the literature yield a finite Griineisen ratio l^ . 

In this Letter we communicate the first-ever observation of a divergent Griineisen ratio 
r for T — ► 0. CeNi 2 Ge2 is known to be a NFL compound which exhibits a paramagnetic 
ground state jl^. The electrical resistivity, p(T), resembles that of CePd2Si2 at the pressure 
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tuned QCP Q: P ~ Po oc T e with 1.2 < e < 1.5 below 4 K 2, IM llM [l6j, ll2| . In YbRh 2 Si 2 
pronounced NFL effects, i.e. C /T oc — log(T) and Ap ~ T have been observed upon cooling 
from 10 K down to 0.3 K. While Ap(T) keeps following the linear T-dependence down to 

n n 

Tjy =70 mK, C/T diverges stronger than logarithmically below T = 0.3 K [5J, |18[ . For our 
study we chose a high-quality single crystal of YbRh 2 (Si .95Ge .o5)2 for which T N has been 
reduced to 20 mK. Large CeNi 2 Ge 2 single crystals of weight 5 to 6 g were grown using the 
traveling-floating- zone mirror-oven technique. The samples were analyzed via electron-probe 
microanalysis (EPMA) and found to have the desired 122 (P4/mmm) structure with little 
mosaic spread, good stoichiometry and no second phases. Two thin oriented bars with a 
length of 5 mm were formed by spark-erosion from the center of the large single crystal. The 
bars were annealed for 120 hours at 800°C in an Ar partial pressure (10 — 20 torr) and have 
a residual resistivity of 3 — 5 /if2cm. Single crystalline platelets of YbRh 2 (Sii_ x Ge x ) 2 with 
a nominal Ge concentration of x = 0.05 were grown from In flux as described earlier 
From a careful EPMA, the effective Ge-concentration is found to be x e // < 0.02 ± 0.01. 
The large difference between nominal and effective Ge-content is due to the fact, that Ge 
dlSs oive S b e tter than Si in tbe Mu , A * effective Ge eonten. of 0.02 ± 0.004 Q 
is deduced from hydrostatic pressure experiments Q|. The residual resistivity of the Ge- 
doped crystal is 5 fiQcm. The thermal expansion and the specific heat have been determined 
in dilution refrigerators by utilizing an ultrahigh resolution capacitive dilatometer and the 
quasi-adiabatic heat pulse technique, respectively. 

Figure 1 displays the T dependence of a a and a c , the linear thermal expansion coefficients 
of CeNi 2 Ge 2 measured along the tetragonal a- and c-axis. As shown by the solid lines, the 
data can be described in the entire T range investigated by the T-dependence predicted jlo| 
by the three-dimensional (3D) SDW scenario, i.e., the sum of (singular) square-root and 
(normal) linear contributions. The corresponding fit parameters are listed in Table I. We 
observe a moderate anisotropy a c ~ 1.8a a . As shown in the inset, the volume expansion 
coefficient (3 = 2a a + a c , plotted as (3(T)/T is not a constant upon cooling, as would be for 
a Fermi liquid, but shows a 1/y/T divergence over more than two decades in temperature 
from 6 K down to at least 50 mK. This is one of the cleanest observations of NFL behavior 
in a thermodynamic property made in any system so far. 

We next consider the low temperature specific heat of CeNi 2 Ge 2 . As shown by several 
investigations, C(T)/T strongly increases upon cooling from 6 K to 0.4 K jj, 
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This increase has either been described by C(T) /T oc — log(T) |4j,lli| or C{T)/T = jq — cVT 
j^j]. Below 0.4 K different behaviors have been reported. While Knopp et al. found a peak 
at 0.3 K followed by a 6% decrease in C(T)/T from the maximum value 



Koerner et al. 

observed a leveling off in C(T)/T below 0.3 K Q|. In contrast, C(T)/T of a high-quality 
sample with very low residual resistivity does not saturate but shows an upturn at the lowest 
temperatures |22j. Very recently a systematic study of the low temperature specific heat 
on different high-quality polycrystals, prepared with a slight variation of the stoichiometry 
jl^ , has been performed. The result was that nearly all of the different investigated samples 
showed such an upturn in C(T)/T below 0.3 K whose size, however, is strongly sample 
dependent even for samples with similar p(T) and a residual resistivity of only 0.2 juQcm 
231 ] . In the following we analyze the specific heat (Fig. 2) measured on the same sample 
that has been used for the thermal expansion study. Below 3 K the data can be described 
by C(T)/T = 70 — cyT + d/T 3 using the parameters listed in Table I (solid lines in Fig. 2). 



Here we assume that the low temperature upturn, present in this sing 



e crystal as well, 



24j . Its influence on 



could be ascribed to the high-temperature tail of a Schottky anomaly 
the Griineisen ratio is smaller than 5% at 0.1 K and therefore not visible in the T(T) plot 
shown in the inset of Fig. 2. This is the first observation of a divergent T(T) for T — ► in 
any material and provides striking evidence that CeNi 2 Ge2 is located very close to a QCP. 
The observed T dependence is in full agreement with the 3D SDW prediction 10]. If the 
investigated high-quality single crystal would enter a Fermi liquid regime below 0.3 K as 
observed for the sample studied in [161 ]. T(T) should saturate below that temperature. 

The application of magnetic fields to CeNi2Ge2 is found to gradually reduce the low-T 
specific heat coefficient. For B>2Ta nearly temperature- inde pen dent 7(5) = C(T, B)/T 
is observed at low temperatures with 7(1?) = 7o — const\J~B |2lJ. The low-temperature 
thermal expansion shows a similar field-induced crossover to Fermi liquid behavior (Fig. 3) 
and the field dependence of a(T, B)/T in the field induced F. 
shown). Both features are consistent with the predictions 
fluctuations at a zero-field AF QCP, assuming a linear dependence between the magnetic 
field and the distance r from the QCP. 

We now turn to YbRh 2 (Sio.95Geo.o5)2, i n which we have measured the thermal expansion 
from 50 mK to 6 K. Compared to CeNi2Ge2, here the volume thermal expansion coefficient 
(3{T) has an opposite sign reflecting the opposite volume dependence of the characteristic 



regime diverges like 1/yB (not 
lp| from the itinerant 3D SDW 
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energies. At T > 1 K, (3{T) can be fit by -Tlog(T /T) with T « 13 K (see Fig. 4). At 
T < 1 K, the best fit is given by a\ + a^T. Both are not only different from the expected 
3D-SDW results discussed earlier, but also weaker than the lnlnT form expected in a 
2D-SDW picture Q- 

The difference from the 2D-SDW picture is even more striking when 
we look at the Griineisen ratio. In Fig. 4, we have also shown the electronic specific heat at 
zero magnetic field. Here C e i = C — Cq, where Cq oc 1/T 2 denotes the nuclear quadrupolar 
contribution determined from recent Mossbauer results |27f. At 20 niK a maximum in 

n 

C e i(T)/T marks the onset of very weak AF order |9|. This is suppressed by a tiny critical 
magnetic field of B c = 0.027 T applied in the easy plane. At B = B c , a power law divergence 
C e i(T)/T oc T -1 / 3 is observed (which is already incompatible with the 2D-SDW picture) 
At higher temperatures the zero-field specific heat coefficient also varies as log(TQ/T) with 
Tq = 30 K (Fig. 4) 5]. Because of the difference between Tq and T , the Griineisen ratio is 
strongly temperature dependent. Below 1 K it diverges as r(T) = r + cT~ 2 / 3 , i.e., weaker 
than the ^ ^t^P form expected in a 2D-SDW picture [l(|. We note that, in the measured 
temperature range, the zero-field data of both the specific heat and thermal expansion are 
identical to their counterparts at the critical magnetic field. 

To interpret our results, we introduce a Griineisen exponent x in terms of the critical 
Griineisen ratio T cr oc ^ oc ^f, where (3 cr and Cy are the thermal expansion and specific 
heat with the background contributions subtracted: this exponent is equal to the dimension 
of the .est relent operator that is coopied to pressure Q. It is show, in Re, Q that, 
for magnetically three-dimensional systems without frustration the SDW picture should 
apply. This is consistent with our finding here that both the thermal expansion and specific 
leat results in CeNi 2 Ge2 can be fit by the respective expressions for a 3D-AF-SDW theory 
3(j] . Our results correspond to f3 cr oc VT and Cy oc T 3//2 , leading to r cr oc h. In other 
words, the Griineisen exponent x = 1 (with error bars +0.05/ — 0.1, as determined from 
a log- log plot shown in the inset of Fig. 3.) In an SDW picture, the most relevant term is 
the quadratic part of the 4 -theory. The corresponding dimension is l/i/z = 1, where the 
spatial-correlation-length exponent v = 1/2 and the dynamic exponent z = 2. 

For YbRh 2 (Sii_ x Ge x .)2, on the other hand, the measured Griineisen exponent is frac- 
tional: x = 0.7 ± 0.1 as determined from a log-log plot of r cr versus temperature shown 
in the inset of Fig. 4. While definitely not compatible with the itinerant SDW theory, a 
fractional Griineisen exponent is consistent with the locally quantum critical point. One 
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kind of condition favorable for this new type of QCP corresponds to a magnetic fluctuation 
spectrum that is strongly anisotropic At such a locally quantum critical point, spatially 
local critical excitations emerge and co-exist with the spatially extended critical spin fluc- 
tuations. There are then two scaling dimensions to be considered. For the tuning of the 
long- wavelength fluctuations, the dimension of interest is still given by the expression 1/vz. 
While v remains 1/2, the dynamic exponent z becomes 2/ a > 2 where a is the fractional 
exponent that characterizes the dynamical spin susceptibility. As a result, 1/vz < 1. For the 
tuning of the local fluctuations, the corresponding dimension is the inverse of the temporal- 
correlation- length exponent. Within an e— expansion scheme as carried out in Ref. |2£| and 
for the XY-spin-invariant case of relevance to YbRl^Sii^Ge^, this exponent is found 2t| 
to be 0.62 to the first order in e and 0.66 to the second order in e. The overall Griineisen 
ratio will then display a fractional exponent, as indeed seen experimentally. 
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Chemischen Industrie (Dresden), the Dutch Foundation FOM-ALMOS (O.T. and J.A.M.), 
NSF, TCSAM, and the Welch Foundation (L.Z. and Q.S.) for support. 
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a(T) = aVT + bT a \\ c a = 1.5 • l(r 6 K- L5 ,& = 0.87 • lO^K" 2 

a || a a = 0.99 • lO^K -1 - 5 ,** = 0.42 • lO^K" 2 
0(T) = aVT + bT (3 a = 3.5 • 10- 6 K- L5 , 6 = 1.7- lO^K" 2 

C(T)/T = 70 - cVf + d/T 3 70 = 0.46JK- 2 mor\c = O.llJmor^- 5 / 2 

d = 102/iJKmor 1 



TABLE I: Fit forms and parameters for CeNi 2 Ge2- 
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FIG. 1: Linear thermal expansions of CeNi2Ge2 vs. temperature at B=0. Inset shows volume 
expansion as (3/T vs T. Solid lines are fits as specified in TABLE I. 
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T(K) 

FIG. 2: Specific heat at B=0 as C/T vs T for CeNi2Ge2. From the raw data (dashed line at low T), 
a contribution C n = a/T 2 with a = 102 /iJK/mol has been subtracted giving the low-T open circles. 
The inset shows the T-dependence of the Griineisen ratio T = V m / 'kt • (3/C where V m and kt are 
the molar volume and isothermal compressibility, respectively. Here, we use kt = 1-15 x 10 -11 
Pa -1 as determined from high-pressure lattice parameter measurements at room temperature j^. 
Solid line is a fit as specified in TABLE I. 
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T(K) 

FIG. 3: Thermal expansion of CeNi2Ge2 along the a-axis as a/T vs T at varying magnetic fields. 
Inset shows the critical Griineisen ratio T cr = V m / • (3 cr /C cr as logr cr vs logT (at B = 0) with 
per = p( T ^ _ bT and C cr = C ( T j _ + d / T 2^ uging the parameters i isted in TABLE I. Solid 

line represents T cr oc 1/T X with x = 1. 
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FIG. 4: Electronic specific heat as C e i/T (left axis) and volume thermal expansion as —j3/T (right 
axis) vs T (on a logarithmic scale) for YbRli2(Sio.95Geo.o5)2 at B=0. Solid lines indicate log(To/T) 
dependences with Tq = 30 K and 13 K for C e i/T and —(5/T, respectively. Dotted line represents 
-0/T = a + oi/T with a = 3.4 • 1(T 6 K~ 2 and a\ = 1.34 • 1(T 6 K -1 . Inset displays the log-log 
plot of T cr (T) with T cr = V m /KT ■ P cr /C cr using k t = 5.3 • 10" 12 Pa" 1 H, P cr = P{T) + a T and 
qct _ Solid and dotted lines represent T cr oc 1/T X with x = 0.7 and x = 1, respectively. 
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